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Abstract 

Magnetic properties of graphene with randomly distributed magnetic defects/vacancies are stud- 
ied in terms of the Kondo Hamiltonian in the mean field approximation. It has been shown that 
graphene with defects undergoes a magnetic phase transition from a paramagnetic to a antiferro- 
magnetic (AFM) phase once the temperature reaches the critical point T/y. The defect straggling 
is taken into account as an assignable cause of multiple nucleation into AFM domains. Since each 
domain is characterized by partial compensating magnetization of the defects associated with dif- 
ferent sublattices, together they reveal a super-paramagnetic behavior in a magnetic field. Theory 
qualitatively describe the experimental data provided the temperature dependence of the AFM 
domain structure. 

PACS numbers: 73.22.Pr, 75.30.Hx, 75.50.Dd, 75.70.Ak 
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Along with unique transport characteristics the magnetic behavior of the graphene-based 
materials attracts much attention in recent studies due to the significant interest in fun- 
damental physics and prospective spintronic applications. Particularly, the possibility of 
the band gap manipulating in antiferromagnetic ordered defective graphene (i.e. graphene 
with adatoms or vacancies)^ - - offers an additional control for nonlinear functionality while 
an efficient spin injection capability into graphene^ with long spin coherence time/length 
even at room temperature^ puts the graphene in the forefront of the materials for emerging 
spin-based information processing. Furthermore, the room temperature weak ferromag- 
netism (FM) has been reported in highly oriented pyrolitic graphite irradiated by proton 
beams& and in defective graphene prepared from soluble functionalized graphene sheets.- 
Moreover, it was experimentally discovered the coexist of ferromagnetic correlations along 
with antiferromagnetic (AFM) interactions in all series of the multilayer defective graphene 
samples in Ref. {]■ On the other hand, no ferromagnetism has been detected in pure graphene 
nanocrystals in wide range of temperatures^ revealing ambiguity about the graphene edge 
contribution to magnetization data (compare Refs. llCH12[). 

It was recently theoretically and experimentally realized that the vacancies or hydro- 
gen addatoms associated with carbon atoms mediate the local magnetic moments in the 
graphene.—"— Moreover, the local spin moments of the defects reveal strong exchange inter- 
action with delocalized electrons that is the main source to trigger low temperature anomaly 
in conductivity and Kondo effect.—"— Theoretical analysis of indirect interaction through 
the graphene carriers in terms of Kondo Hamiltonian unambiguously indicates the ferromag- 
netic (antiferromagnetic) exchange coupling between localized spins associated with the same 
(different) sublattices of the crystalline graphene.—"— This fact as well as Lieb's theorem^! 
are usually adduced to back up the arguments that a weak ferromagnetism is a result of 
imbalance ANd in the number of vacancies/impurities of the A or B sublattices (i.e. not 
precise compensation of sublattices magnetizations with opposite directions).— 1 ^ At the 
same time the actual imbalance of randomly distributed spin moments fails even a weak 
ferromagnetism because |AA^| /Na — > in large graphene crystal with the total number 
of the defects. 

In present study we propose a different approach to the problem of coexistence of FM 
and AFM phase in a single graphene layer. It takes into account a multiple nucleation of 
antiferromagnetic domain germs in the mass of randomly distributed spins in monolayer 



2 



graphene. At the same time, the domains each reveal a small random imbalance And 
and magnetism so that together they additively contribute to the net magnetic moment 
culminating the saturated magnetization in a magnetic field. In such a way the finiteness of 
the magnetic correlations should be taken into account. The strong short-range correlation 
within each domain is responsible for its antiferromagnetic ordering. Such AFM correlations 
are weakened or broken at the domain boundaries, which can be associated with the impurity 
rarefactions or lattice imperfects. Thus the behavior of the magnetic domains in a magnetic 
field is almost mutually independent that constitutes some finite magnetization in the limit 
— > oo. It is important to note that magnetism in graphene does not come into conflict 
with Mermin- Wagner theorem^ because of finite magnetic anisotropy^ 5 ^ and restricted 
sizes of AFM domains. 

The microscopic approach is based on the Kondo Hamiltonian of carrier-localized spin 
exchange interaction that has been extensively used in the analysis of magnetic properties of 



graphene in Refs 



.y 
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21 



and 



26l . The spectrum of bulk graphene is treated in tight binding 
approximation. In mean field approximation, this approach leads to set of two equations 
for the sublattice magnetizations that define a critical temperature of the antiferromagnetic 
ordering without any limitations on the number of carbon atoms involving to computation 
in more sophisticated models (Refs. Q and 0). As a result, we find a simple analytical 
expression for Neel temperature that may serve as a guide for analysis of the numerous 
experimental data on graphene magnetism and estimate the mean size of the magnetic 
domains in graphene. 

Let us consider a graphene fragment possessed large enough a flat area Af to neglect the 
edge effects. The Hamiltonian in momentum representation takes the forrn^ 

#k = 7ccpi/i(k) + a 2 / 2 (k)], (1) 

where 7 CC = 2.7 eV is the matrix element of electron hopping between nearest neighbor 
atoms connected with the vectors e m (m— 1,2,3), the Pauli matrixes cXj are defined over the 
sublattices A and B basic functions and k =(k x ,k y ) is the electron momentum; /i(k) = 
X] m cos(ke m ), / 2 (k) = X!m sin ( ke ™)- We assign e m = a cc (cosmw,sinmw), where a cc = 
0.142 nm and u = 2ir/3. In diagonal form, Hamiltonian (TjQ) describes the graphene dispersion 
law 

e&,k = lccbe(a k) (2) 



for conduction (6 = 1) and valence (b = — 1) bands with e(aok) = a/ ff (k) + /|(k). A 
straightforward algebra shows that e(q) = ( 3 + 4 cos cos y + 2 cos q y ) , i.e. there 
are two non-equivalent contact points at the corners of first Brillouin zone (BZ) K, K' = 
^■(a/3, ±1,) where conduction and valence bands are degenerated at e^K = %k' = 0; 
a = \/3a cc is a length of lattice vectors. In the vicinity of these points the energy dispersions 
are the Dirac cones, e^ K (,) +)( = bhvp with Fermi velocity Vp = -y7 cc a o/^- 

The Kondo Hamiltonian of the exchange interaction between a band electron (with po- 
sition r and spin S) and the n d localized spin moments lj pinned to the sites Rj of the 
graphene lattice reads 

H K = -J2 J ^ R ^ ( 3 ) 

where J(r, Rj) rj Jwq5(t — Rj), J is the exchange constant, = ua + tib the total number 
of the defects located at A and B sublattices, w Q = \f?>a^/2. In the bipartite graphene 
lattice, it is convenient to double group the Hamiltonian of Eq. ([3]) on sublattice defects 
3a{b) — I; •■■> n A{B)'- Hk = + In the representation of eigenfunctions of Eq. (Op), 

each part of Hk manifests itself through the projection operators Pa(b) — (1 ± <?3)/2. Then 
the summing up over the random scattered impurities and thermal averaging of their spin 
states reduce Eq. (j3J) to 

H ex = a(mS + wC), (4) 

where the total magnetic moment of the defects m = + mg and their antiferro mag- 
netic vector w = - tng are expressed via sublattice magnetic moments m^^) = 



—uaQ^b (li A (B))j a = nj~g/l^> n f = Af/vjQ the total number of primitive sells, zu is their 
area, g = 2 and \xb are the <?- factor and Bohr magneton. The distinguishing property of 
Eq. (pi]) is that AFM vector exerts the composite spin C = a$ <g) S that can be a finite 
magnitude even under (S) = 0. 

Let us introduce the effective magnetic fields and B b that cause the spin polarizations 
in each sublattice (Ij A(B) \ = ~^a{b) f tanh 9IJlB ^ {B) , n A (B) = &a(b)/Ba(b), k B = 1. On the 
other hand they can be represented as 33 

B.(B)=-^- (5) 
in terms of of thermodynamic potential 

$ = -T ^ In (1 + e^~ E ^ T ) , (6) 



and chemical potential /x. The energy bands of the defective graphene with Hamiltonian 
Hk + H ex are described by 
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E b , s ,u = ^am + b X /El k +(— J , (7) 

where s = ±1 is a spin number. A distinctive feature of the dispersion law ([7]) consists 
in opening bandgap E g = a(w — m) provided that w > m. Such splitting of the energy 
bands lowers the total electronic energy of the valence band that cannot be compensated 
by raising electron energy in the conduction band resembling the cooperative Jahn- Teller 
effect. 34 Equation ([5]) constitutes the closed set of the equations for both = zBa and 
Bg = zBb, the unit vector z is directed along quantization axis. Explicitly they take the 
form of integral over the first BZ, 



&A,B — 



[ [ dq x dq y F(k); (8) 

J J BZ 

*w = j > : i o ± , baw n i m >s , k - a*), (9) 

4 A 2 k +(T)' 




where /(i£&,s,k — fi) is the Fermi-Dirac function and ± discriminates the different sublattices. 

In the first stage we focus on the stronger effect of AFM ordering within a single domain 
that eventually establishes weak FM. Ignoring small imbalance And <C rid, Eqs. ([8]) and (Q 
decompose on independent equations for Ba and Bb = —Ba so that each Ba and Bb satisfy 
to the self-consistent equation 



J 2 x 



J(y,/i,r)taohf^, (10) 



in terms of variable y = gfisBA(B)/T. A complex integral function I(y, /x, T) can be reduced 
with high accuracy to the constant Iq ~ 0.448 provided that 7 CC ^> Jx,fi,T. Therewith in 
the limit y — > Eq. (flOj) gives rise the expression for Neel temperature 

J 2 x 

T N = ^I , (11) 

which depends merely on the exchange constant and defect fraction x = nd/2rif. The squared 
dependence of Tjy on J as well as the quantitative evaluation T/v = 0.01127 cc at x — 0.2 
and J = 0.57 cc are in a good agreement with Monte-Carlo simulations carried out in Ref. 
I2I. Particulary, the Eq.t flT]) shows that exchange constant |J| = 1.9 eV guarantees a room 
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temperature AFM ordering at reasonable low x = 0.04. Solutions y = y(x,T) of Eq. (flQj) 
define also the energy gap E g = Ja; tanh y ^ T ^ and magnetization tua = \riAg^B tanh y ( x ^ , 
Fig. 1 illustrate the temperature dependence of E g calculated with Eq. ( IT01) . 

Nucleation of AFM domains each in possess of a random imbalance ha — ns lays 
the groundwork for weak ferromagnetism in graphene. Applying the statistical approach, 
let us attribute x to the equal probability for each site A or B to be defected for all 



domains so that mean number of deviation is {tia — ub) = because tia = tib = xrif. The 



mean-square estimate is {tia — tib) 2 = 2njx(l — x) = nd(l — x), thus the mean number of 



non-compensated spins is Ana = ^/2n/x(l — x) = \fndil — x). The respective magnetic 
moment 



m(nf, x, T) ~ fiB y 2nfx(l — x) tanh — ^ — (12) 

can be interpreted as a weak ferromagnetism attributed to a single domain. 

Above we considered the formation of the magnetic moments m(rif, x, T) < fiBy/^d { n d ^ 
y/n2 ^> 1) in individual AFM domains. To consider the behavior of multi-domain graphene 
in a magnetic field B further conjectures must be done. It is convenient to split the ensemble 
of domains on subsets, each determined by certain space Af and number of the defects. 
In a magnetic field B the Af -domain with defects contributes to the net magnetization 
as m(nf,x,T)L[Bom(nf,x,T)/T e ff], where L(x) is a Langevin function and T e ff = T + Taf 
is an effective temperature. Parameter Taf takes into account the tendency to merge several 
small domains into one large AFM one that reveals a close analogy of this parameter with 
the temperature shift T in diluted magnetic semiconductors with AFM interaction between 
localized spin moments. 35 Such interdomain interaction constitutes proportionality of Taf 
to the length of domain boundary network or inverse proportionality to domain mean size 
L. The final magnetization output M can be expressed in terms of distribution function 
f(Af, x) so that 



M = J J m(n f , x, T)L ^ B o m (n f ,x,T) ^ ^ 



Apparently the dispersion in the domain shapes and spaces and in the defect densities is 
very specific for particular sample preparation, chemical and thermal treatment; therefore 
the f(Af,x) cannot be specified a priori. So in the rest part of the paper we focus on 
the analysis of particular experimental data on graphene magnetization reported in Ref. 
[7. There was found that at room temperature, a magnetic field B ~ 3 kOe saturates 



magnetization at small amounts M s = 0.02 emu/g and M s = 0.004 emu/g in two different 
graphene samples Gr600 and Gr400 respectively. At the same time a stronger magnetic 
field B Rj 30 kOe needs to saturate magnetization with much higher M s = 0.8 emu/g and 
M s = 0.2 emu/g at T = 2 K. Besides, in very narrow magnetic field region there was recorded 
some hysteresis loop stipulated by a weak spin-orbital coupling. This effect is unobtrusive 
in the scale of magnetization curves recorded in Ref. Q and is not discussed here. 

Two different scenarios can be applied to these results. First, let us assume that the 
f(Af,x) is unchangeable in all range of the temperatures. In such a case the experiment 
implies that the majority of the domains possess relatively small area with slight magnetic 
moments m(nf, x, T) and relatively low susceptibility at T = 2 K. With temperature increase 
this portion of the domains undergos to phase transition to paramagnetic state and drop out 
from the consideration. The rest part of large domains with higher magnetic susceptibility 
can be responsible for magnetization at T = 300 K. This scenario, however, fails to describe 
the temperature variations of the curves M = M(Bq) predicting M s to be much smaller 
than the observable magnitudes at T = 300 K. 

Another approach assumes that the parameters n/ = Aj/wq and rid = 2xnf obey 
to normal (Gaussian) distributions G{nf) and G{rid) around their mean values nj and 
rid but nj can vary with T. Substituting f(Af,x)dAfdx = f(nf,nd)dnfdrid, where 
f(nf,n d ) = N D G{nd)G{rif) into Eq. f|T3|) allows to describe all experimental data with 
good accuracy (Fig. 2). At the same time one must assume a growth of the domain mean 
sizes L = Aj with temperature. Note that the related effect of Taf decrease as L _1 (as was 
discussed above) correlates with this model (Fig. 2). Such variations as well as detail domain 
structure are caused by structure inhomogeneous that is beyond the developed theory. At 
the same time, the strengthening of AFM correlations through the domain boundaries with 
temperature increase seems not surprising if high temperatures favors to electron overcome 
the barriers between different domains. Note also that the intriguing result of vanishing 
ferro magnetism in the sample^ Gr800 might just be the effect of inhomogeneous removing 

— 1/2 

after high temperature annealing that decreases magnetization as Af — > 0. 

In conclusion, we have shown that the imbalance in the numbers of the defects located 
at A and B sublattices along with graphene fragmentation into the AFM domains results in 
small but finite magnetization response in a magnetic field. Theory quantitatively describe 
experimental data provided that the domain structure can vary with temperature. The 
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following experimental verification needs to establish the model applicability for particular 
samples of the defective graphene. 
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FIG. 1: Band gap of the defective graphene vs temperature at J = 1.9 eV and defect molar 
fractions 4%, 5% and 6% (curves 1, 2 and 3 respectively). 
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Magnetic field (kOe) 



FIG. 2: Magnetization curves of the defective jgraphene calculated with Eq. (|13p (solid and dashed 
lines). Experimental data are taken from Ref. 7 for the samples Gr600 (circles) and Gr400 (squares) 
for T=2 K (a) and T=300 K (b). Fitting to the Gr600 data provides the x = 0.19 and mean sizes 
of the domains L = 43 nm (T = 2 K, T AF = 58 K) and L = 1.6 /mi (T = 300 K, T AF = 1.5 K); 
respectively x = 0.12, L = 130 nm (T = 2 K, T AF = 18 K) and L = 5.6 fj,m (T = 300 K, T A f = 0.5 
K) have been applied to Gr400. 
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